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ABSTRACT 

We use CO and Ha velocity fields to study the gas kinematics in the spiral 
arms and interarms of M51 (NGC 5194), and fit the 2D velocity field to estimate 
the radial and tangential velocity components as a function of spiral phase (arm 
distance). We find large radial and tangential streaming velocities, which are 
quahtatively consistent with the predictions of density wave theory and support 
the existence of shocks. The streaming motions are complex, varying significantly 
across the galaxy as well as along and between arms. Aberrations in the velocity 
field indicate that the disk is not coplanar, perhaps as far in as 20" (800 pc) from 
the center. Velocity profile fits from CO and Ha are typically similar, suggesting 
that most of the Ha emission originates from regions of recent star formation. 
We also explore vortensity and mass conservation conditions. Vortensity conser- 
vation, which does not require a steady state, is empirically verified. The velocity 
and density profiles show large and varying mass fluxes, which are inconsistent 
with a steady flow for a single dominant global spiral mode. We thus conclude 
that the spiral arms cannot be in a quasi-steady state in any rotating frame, 
and/or that out of plane motions may be significant. 

Subject headings: galaxies: individual (M51) - galaxies: spiral structure - galax- 
ies: kinematics and dynamics 



1. Introduction 



Spiral arms are the dominant morphological features of most disk galaxies. From a 
theoretical perspective, two frameworks have been proposed to describe the nature of the 
spiral arms: one is that the spiral arms are generall y long-lasting, or slowly e volving, and the 



other is that the arms are transient features (e.g. iToomre &: Toomrelll972l ). Observational 



studies have yet to show definitively whether the arms are evolving or long lived, though it 
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has been over 40 years since the landmark paper by lLin fc Shul (119641 ) suggesting that spiral 
structure in galaxie s is a long lived phenomenon — the Quasi-Stationary Spiral Structure 



(QSSS) hypothesis (lLindbladlll963l ). In the QSSS depiction, though material passes in and 



out of the arms, the slowly evolving global pattern rotates with a single angular speed that 



Bertin et al. 


1989aljbl: 


Bertin & Lin 


1996) 



However, interaction between a disk galaxy and a companion is another explanation for 
the presence of spiral arms. In such a framework, the arms are transient features that are 



generated by the tidal interaction (e.g.. iToomre &: Toomrdll972l). Any spiral arms e xisting 



before the encounter are overwhelmed by the tidal driving (jSalo fc Laurikainenll2000l ) 



Regardless of the origin of the stellar arms, gas in the disk will respond strongly to 
the gravitational perturbations those arms impose. Numerical studies have indicated that 
shocks can develop if the relative speed between the spiral perturbation and the gas is 
large (iRobertd 1 19691 : IShu et al.l Il973l : IWoodwardI 119751 ) . The presence of dust lanes in the 
spiral arms and the enhancement in ionized emission downstream, indicating regions of 
star formation, is attributed to this shock scenario. Such shocks are also thought to be 
the cause of the well defined molecular arms seen in many grand design galaxies, including 
M51. Numerical and analytical studies have provided predictions for the velocity and density 
profiles of the matter affected by the spiral gra vitational perturbation (e.g 
Kim fc Ostri"k^l2002l : Icittins fc Clarkell2004h . 



Lubow et al. 1986 



There have been numerous observational studies addres sing the nature of spiral structure 
that have focused on the gaseous components. IVissen (119801) showed that ste ady state density 
wave models fit the HI kinematics of M81 quite well. iLowe et al.l (119941) used the m odal 
theory of density waves to describe the spiral pattern in M81. Both iRandl (119931 ) and 
Aalto et al.l (119991 ) used observed molecular velocities along ID cuts on the major and minor 
axes of th e grand design spiral M51, a nd found qualita t ive ag reement with the density wave 



models of [Roberts &: StewartI (119871 ). iKuno fc Nakail (119971 ) fitted observed CO velocities 



from single dish observations to obtain gas streamlines. The smooth shape of the velocity 
pr ofiles led the r a to c onclude that galactic shocks do not exist in M51. However, the study 
by lAalto et al.l (Il999l ). using higher resolution interferometric data, found steeper velocity 
gradients, supporting the presence of shocks. 

Yet, other observational studies have sugg ested t hat t he arms are not long lived. In 
fact, the classic kinematic study of M51, that of iTullyl (119741 ). found evidence for a transient 
pattern in the outer disk, due to the interac tion with its companiori , but that a steady state 
is probably appropriate for the inner arms. lElmegreen et al.l (19891) and IVogel et al.l (119931 ) 
also suggested the presence of multiple pattern speeds. iHenry et al.l (120031 ) argued that the 
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spiral pattern may be a superposition of an m=2 mode and a weaker m=3 mode, suggesting 
a transient pattern for the arms of M51. 

This paper presents a detailed study of the gaseous velocity structure associated with 
the spiral pattern in M51. In a future paper, we will discuss and compare the spiral pattern 
in different tracers. Here, we use the CO and Ha velocities to map the 2D velocity field in 
M51. 

Our study makes use of the full 2D velocity field in M51 from interferometric CO and 
Fabry-Perot Ha observations, rather than just major and minor axis cuts. Noting that 
variations in the observed velocity field are mainly associated with the spiral arms, we fit the 
observed velocity field to obtain the radial and tangential components as a function of arm 
phase (i.e. distance perpendicular to the arm). We then analyze whether the fitted velocity 
field and density maps are consistent with the predictions of steady state theory. 

In the next section we briefly describe our CO and Ha observations. In ^we describe the 
method we employ to estimate the radial and tangential velocity components throughout 
the disk. Since our method is sensitive to the assumed values of the systematic velocity, 
major axis position angle, and disk inclination with respect to the sky, in §1] we present 
results from our effort to constrain these parameters, and describe how errors could affect 
the fitting results. In ^ we present and discuss the fitted profiles of radial and tangential 
velocities for a range of radii. We then use the velocity and density profiles to empirically 
test conservation of vortensity in §6.11 Next, in §6.21 we examine whether (quasi) steady 
state mass conservation is applicable, as would be necessary for a QSSS description. Finally, 
in ^ we summarize our conclusions. 



Observations 



The CO and Ha intensity and velocity maps are shown in Figures [T] and [H respectively. 
The CO J=l-0 data for M51 were obtained in part from BIMA S ONG (Survey of Ne arby 
Galaxies). The ob servations and data reduction are described in iRegan et al.l (120011 ) and 
Heifer et al.l (120031 ). The SONG map is based on 26 pointings and has an angular resolution 
of 5.8" X 5.1". Later, we obtained data for 34 additional pointings, so that the spiral arms 
were mapped as far as the companion galaxy to the north and to a similar distance along 
the spiral arm to the south. Additionally, inner fields were mapped in a higher resolution 
array (B array), yielding higher a ngular resolu t ion. T he newer data were reduced using the 
same procedures as described in iHelfer et al.l (120031 ) for BIMA SONG. Together the data 
sets cover 60 pointings. The maps used for this paper have variable resolution, reaching as 
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high as 4" in the inner spiral arms but degrading to 6" - 13" in the interarms and in the 
outer arms. 

Ha data were obtained with the Maryland- Cahech Palomar Fabry-Perot, which covered 
the optical disk at an angular resolution of 2" and a velocity resolution of 25 km s~^. The 



Gruendl 


( 


1996) and also 


Voeel et al. 


( 


1993) 



Both CO and Ha intensity maps are obtained by fitting Gaussian profiles to the spectrum 
at each location. The velocity maps indicate the velocity of the peak of the fit Gaussian 
intensity. 

Also shown in Figures [1] and [2] are two lines tracing logarithmic spirals. The bright CO 
arm is well represented by a logarithmic spiral with a pitch angle of 21.1°. The weaker arm 
also generally follows a logarithmic spiral, although, as will be discussed, a number of its 
arm segments either lead or lag the depicted line. The logarithmic spirals will be discussed 
extensively in the following sections. 



3. Estimation of Spiral Streaming Velocities 



The observed line of sight velocity Vobs can be decomposed as a sum of terms involving 
the systematic velocity Vsys, the radial velocity vr and the tangential velocity vg: 



Vobs{R, 9) = Vsys + [vr{R, 0) sin(^ - Oma) + ve{R, 9) cos{9 - 9m a)] sini, 



(1) 



where R and 9 are the galactocentric radius and azimuthal angle, and 9ma and i are the 
position angle of the major axis and inclination of the galaxy, respectively. This equation 
does not include a velocity component perpendicular to the disk. The exclusion of the vertical 
velocity component is reasonable since studies of face-on grand design spirals ind icate that 
the 2;-component of velocity is less than 5 km s~^ flvan der Kruit fc Shostaklll982l ). provided 
the disk has no significant warp (we return to this issue in §4.21) . 

Inspection of the velocity maps indicates that the isovelocity contours near the spiral 
arms tend to run parallel to the arms. For a disk in pure circular rotation and with a fiat 
rotation curve, on the other hand, the isovelocity contours of the projected velocity field are 
purely radial. It is evident that the velocity field of M51 is significantly different from this 
sort of simple "spider diagram," due to the non-axisymmetric perturbations associated with 
spiral streaming. Clearly, vr and Vg vary with azimuth. 

Previous estimates of streaming velocities have used observed velocities near the major 
axis (where the projections of vr vanishes) to estimate vp , and veloci ties near the minor axis 
(where the projections of vg vanishes) to estimate vr (e.g.. lRandlll993l ). However, much of the 
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CO gas is organized i nto G MCs and larger complexes known as GMAs (IVogel et al.lll988l ). 
Further, lAalto et al.l (119991 ) have found that the streaming velocities of M51 GMAs in the 
same spiral arm have a significant dispersion. Observations along a single cut, e.g. the major 
axis, sample discrete GMCs and therefore may give a misleading estimate of the streaming 
velocities. As an alternative, an approach that fits a streaming profile to all the observed 
velocities in an annulus as a function of distance from the arm peak may better characterize 
the streaming velocities. Also, the gas surface density distribution varies significantly at any 
distance from the peak (i.e. the gas is clumpy), and so averaging parallel to the arm may 
better characterize the variation in gas surface density as a function of arm distance. 

Typically, 2D fits to a galaxy velocity fie ld assume that f b and vq are constant along 
rings (e.g. the tilted ring analysis described in iBegemanl (119891 )). By contrast, as mentioned 
earher, vr and vg do vary with azimuth, and indeed inspection of Figures [H and [2] indicates 
that the primary variations are due to the flow through the spiral arms rather than variations 
with galactocentric radius. As in most galaxies, the rotation curve of M51 is relatively flat; 
the radial variations that do occur are associated with spiral arm streaming. Thus, we are 
motivated to assume that radial variations of azimuthally-averaged quantities are negligible 
(at least over relatively limited radial ranges) and that vr and vq vary primarily with spiral 
arm phase ifj. The left panel of Figure [3] shows the relevant geometry depicting the spiral 
arm phase. Our assumption is that vr and vg are constant along narrow spiral arcs, such as 
the segments in Figure [HI that are congruent to the spiral arms. Thus, we rewrite equation 
([T]) (for a limited range of radii) as 



Vobs = Vsys + [vR{tp) sm{9 - Oma) + ^e(^) cos(6' - 9ma)] sini 



(2) 



In order to simplify the process of identifying regions of constant arm phase, we adopt 
a coordinate system in which the spiral arms are straight. lElmegreen et al.l (119891 ) show 
that the spiral arms of M51 appear as straight line segments in a [9, log{R)), or logarithmic 
polar, coordinate system. The right panel of Figure [3] shows the logarithmic polar diagram 
corresponding to the features in the left panel. Figure H] shows the CO intensity and velocity 
maps of M51 in log-polar coordinates, and Figure shows the corresponding Ha maps. 

The sky images in Figures [1] and [2] are first deprojected before being transformed into 
a {9,log{R)) coordinate system. In order to deproject the sky view of a galaxy, the center 
position, position angle, and inclination are required. We initially use the canonical values for 
these parameters, which are listed in Table [H We discuss the estimation of these parameters 
in the next section. 



The two straight lines overlaid on Figures H] and [5] indicate the adopted pitch angle of 
21.1° and also correspond to the spiral loci shown in Figure [H It is clear from the overlaid 
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lines, which are sepa rated by 180° , that the weaker arm is not symmetric with the brighter 



one, as discussed by iHenry et aLl (120031 ) . Yet, both CO arms wrap around approximately 
360° of the galaxy, even though they appear to jump in phase at one or more positions. 
The arms in Ha show more jumps in phase and variations in the pitch angle. In spite of 
their asymmetries, the CO spiral arms are particula rly well described as log arithmic spirals. 



better even than Ha, or the optical arms shown by lElmegreen et al.l (119891 ). 



We will refer to overlaid logarithmic spiral arcs (or lines) as "slits," for we will extract 
observed CO and Ha velocities as a function of position along the arc (or line), similar to 
obtaining long-slit spectra. Each slit marks a region of constant arm phase ip. Thus, while 
the observed velocity varies along the slit due to projection, vr and ve are assumed constant. 
We arbitrarily define the arm phase marked by the leftmost slit in Figures H] and as = 0°. 
The other CO arm appears at an arm phase of approximately ip = 180°; other features such 
as the stellar arms and the gravitational potential minimum may of course be offset from 
the CO arms. 

As noted previously, our fit will assume that the intrinsic vr and vq are constant at 
a given arm phase, i.e. along a given slit, but that vr and vg vary with ip as the slit is 
translated in azimuth. Translating the slit amounts to shifting a straight line to the right 
in the logarithmic polar diagram; this direction of increasing azimuth is the same as the 
direction of rotation for M51. We then fit equation ([2]) to the observed velocities extracted 
at each arm phase ip, thereby obtaining vr{iI)) and vg{tp). 

Although Vr and ve vary primarily with arm phase, they may of course also vary with 
radius. Therefore, we limit the radial range of an annulus (or equivalently the length of a 
slit) as much as possible while still fitting a sufficiently extended azimuth range to obtain 
good leverage on both vr and vg. In other words, an annulus should be sufficiently broad to 
cover the spiral arm both near the major axis and the minor axis; the width of the annulus 
thus depends on the pitch angle and galactocentric radius of the arm. 

We first test our method by applying it to a model spiral galaxy with known radial 
and tangential velocities. The solid lines in Figure [HI shows the averaged density, vr, and vg 
profiles in an annulus from a snapshot of a hydrodynamical simulation of a disk responding 
to a spir al perturbation . The model spiral galaxy is a 2D version of a 3D model described in 



detail in lGomez fc Coxl (120021 ). and the annulus used here extends from 8.38 - 8.92 kpc. The 
direction of gas fiow is in the direction of increasing phase. As the gas approaches the arm, 
(marked by density maxima) the radial velocity, vr, decreases by ~40 km s~^. The sign 
reversal of vr indicates that the gas is moving away from the nucleus before the shock and 
towards the nucleus after the shock. As the gas emerges from the arm, the radial velocity 
increases again. The tangential velocity vg gradually decreases as the gas approaches the 
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arm, then receives a strong boost and reaches a maximum just downstream from the arm. 

In order to test the fitting algorithm, the model vr and vq at all locations are used 
in equation ([T]) to create a model observed velocity field. This velocity field, along with 
the model density map, are transformed into logarithmic polar projections. Equation ([2]) 
is then fit to the model observed velocities at each arm phase in an annulus, using slits 
parallel to the spiral arms; the dashed lines in Figure [6] are the results of the velocity fits, in 
the same annulus (8.38 - 8.92 kpc). The results reproduce the overall shape of the velocity 
profiles quite well, although with slight phase shifts and offsets. The offsets and the shallower 
minimum in vr are likely due to the variation of the pitch angle of the arms with radius; i.e. 
the spiral arms are not perfectly logarithmic, whereas the "slit" used to extract velocities at 
constant arm phase is. Despite these offsets, we were able to reproduce the major features of 
the velocity profiles of the model spiral galaxy, indicating that our method of fitting observed 
velocities at each arm phase recovers a 2D velocity field reasonably accurately. 

We now apply our fitting method to the M51 data, adopting systematic parameters 
listed in Table [TJ As an example. Figure [7] shows the vr (top) and vg (bottom) fits to the 
observed CO and Ha velocity field for one annulus between the galactocentric radii of 21" 
and 36". The CO intensity averaged along a slit as a function of phase angle ip is also 
shown as dashed lines, indicating the distribution of molecular gas. As mentioned, ip = 0° is 
arbitrary and is marked by the leftmost line in Figures H] and O corresponding to the brighter 
arm, which we will refer to as Arm 1. We show a phase range greater than 360° so that both 
upstream and downstream velocities can easily be seen for both arms. The direction of gas 
fiow through the arms (assuming we are inside corotation) is from left to right, so that the 
right sides of the CO peaks correspond to the downstream side of the arm. In most cases, 
as the gas fiows through the arm, the radial velocity decreases and then increases, and the 
tangential velocity receives a boost, as predicted qualitatively by density wave theory. 

In a conventional tilted-ring velocity fitting analysis, galaxy parameters such as the 
inclination, position angle, dynamical center, and systematic velocity can be directly fit. 
However, even though the inclination and position angle appear explicitly in equation ([2]), 
for our fits all but the systematic velocity must be assumed prior to deprojecting a galaxy 
velocity field and therefore before the fit. For our initial fits we employed the standard 
assumed values for these parameters for M51, shown in Table [H In the next section we 
explore the effects of errors in these assumed global parameters on the estimation of vr and 
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4. Method Testing and Parameter Constraints 

To test the sensitivity of the vr and ve fits to errors in the global (fixed) parameters, we 
generated test velocity fields and created sky projections with known parameters. We then 
applied our fitting technique to estimate vr and vg for the model galaxies, assuming incorrect 
values for the fixed parameters, and compared the fitting results to the actual model values 
of Vr and vg. This enables us to quantify the sensitivity of the fits to the fixed parameters. 
In addition to constraints obtained from fitting our kinematic data, we also use standard 
methods to constrain the values of Vgys, Oma, and i. As we shall show, one of our conclusions 



is that some of the basic parameters for M51, many of which date to iTuUyi (119741 ). may in 
fact be poorly constrained due to the morphological and kinematic perturbations induced 
by the tidal interaction with its companion. 

As an initial test, we generated a simple model with vg = 240 km s~^ and vr = —35, 
i.e. an axisymmetric disk with a fiat rotation curve and uniform radial infiow. We refer to 
this model as the "constant velocity" model. We apply our general method to fit vr and vg 
as a function of arm phase using "observed" velocities. If we assume the input values of Vgyg 
= 464 km s~^, 9m A = 1T0°, and i = 20°, we indeed recover the input values for vr and vg 
as independent of phase. We now consider the effects of assuming incorrect values for the 
parameters. 



4.1. Position of Dynamical Center 



In testing the sensitivity of the fits to the assumed center position, we applied the fitting 
algorithm to a model for which the center position was shifted by 1" in both RA and DEC. 
We found that a 1" error in the assumed center has a negligible effect on the fit velocities. 
BIMA observations have an astrometric accuracy of ~10% of the synthesized beam. The 
highest resolution of our CO observations is 4.5", so the error in position will likely not be 
greater than ~0.5". Thus, observational errors will likely not affect the results of our fits. 
For all the analysis that follows, we will adopt the center position listed in Table [TJ This 
choice assumes that the dynam ical center coin c ides with the location of a weak AGN known 
to exist in the nucleus of M51 (IHo et al.l 119871 : iNakai &: Kasugalll988l ). We use the position 
of the radio continuu m source observed with the VLA, which has an accuracy of itO'.'Ol 
( IHagiwara et al.ll200ll ). 
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4.2. Systematic Velocity 

Before discussing mettiods for determining Vgys, we first explore the effect that an error 
in an assumed value of Vgys would have on fits for vr and ve in which Vsys is held fixed, using 
the constant velocity model. As expected, an error AVgys in the assumed Vgys produces a 
sinusoidal variation in both fitted velocity components, with an amplitude of AVgys/ sini, 
and a period of 360° (see eqn. [1]). Clearly, Vsys needs to be well determined. 

One approach to obtaining Vsys is to fit the data for its value using equation ([2]). Figure 
E] shows the results of fits to the M51 data in which Vsys was fit, along with vr and vg, as a 
function of arm phase. Although Vsys should be constant, it can be seen that the fit value of 
Vsys varies with arm phase. Similar variations in fitted Vsys result regardless of what values 
of the position angle and inclination are assumed. 

One possible explanation for the apparent variation of Vsys is that the galactic disk 
of M51 is twisted and/or warped, i.e. the position angle and/or inclination may vary with 
radius. We therefore compare the results of fitting for Vsys from a model galaxy with no warp 
to a model with a warp. We again make use of the constant velocity model, for this model 
also represents an unwarped disk. Instead of keeping Vsys fixed, we allowed this parameter 
to be free in the fit. If we use the true position angle and inclination, we correctly recover 
the adopted values for all three free parameters, Vsys, vr, and vg. 

To generate a warp model, we increase i monotonically from 25° at the inner radius 
(100") to 35° at the outer radius (200"); the position angle is kept fixed. If we allow Vsys 
to be free in fitting the model data. Figure M shows that the Vsys varies almost sinusoidally 
about the true model value 400 km s~^. The mean fitted systematic velocity is equal to 
this value. This is the case regardless of what values of inchnation or position angle we use 
(within reasonable limits), and regardless of the radius range of the annulus used for fitting. 
Thus, regardless of the assumed fixed parameters and of the limits in radius, for a simple 
warp the mean value of the fits gives the correct systematic velocity. 

Motivated by our finding that even with a warp the average fit value of Vsys gives the 
true value, we calculated the mean of the Vsys values shown in Figure [8], obtaining Vsys = 
470.6 km s~^. Comparison of the Vsys fits of the actual M51 data (Figure [8]) with Vsys fits to 
the simple warp model (Figure [H]) shows that the warp model has a slower variation. Hence, 
if a warp is responsible for producing variations in the fitted systematic velocity it must be 
more complex than our simple model; we return to this question in §6.2[ 



We therefore apply two additional methods to estimate the v alue of Vsys- The first 
method is based on a standard tilted-ring analysis (jBegemaru 119891 ) . in which the galactic 
disk is represented as a series of nested tilted rings. In its most general form, each tilted ring 
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may have a different center, systematic velocity, position angle, inclination, and rotational 
velocity. We use 10" rings from an inner radius of 20" to an outer radius of 120", fixing the 
center position, inclination, and position angle to the values in Table [H We obtain mean 
systematic velocity of 471.4 ± 0.5 km s~^. When we allow the position angle to vary as well, 
we obtain a mean of 471.3 ± 0.3 km s~^. 

A third method we use to constrain V^^^ is to assume a functional form for the rotation 



curve, using the NEMO program rotcurshape (jTeubenlll995l ). In contrast to the tilted-ring 
method, which fits each ring independently, rotcurshape fits Vsys, Oma, h center position a 
and 5, and the coefficients of the function used to describe the rotation curve simultaneously 
to the entire velocity field. Therefore, it can yield a single Vgys that best fits the entire 
velocity field. It is particularly useful for finding Vsys if the kinematic center position can 
be fixed. For this fit, we limit the rotcurshape fit to the inner 20" in radius. This is inside 
the main spiral arms, in the region where the rotation curve is rising and the isovelocity 
contours are relatively straight. We assume the center, 6 ma, and i listed in Table [1] and a 
rotation curve of the form v = Vsys + Vox/{l + x) where x is the ratio of the radius to the 
core radius, and fit for Vsys, Vo, and the core radius. We obtain Vsys = 473 ± 0.5 km s~^. If 
we allow the position angle to vary as well, we obtain 473.2 ± 0.3 km s~^. It is encouraging 
that this is within 2 km s~^of Vgys determined from the other two methods even though this 
fit uses a different method and fits an entirely different region (i.e. the inner 20", inside the 
main CO arms, as opposed to outside 20"). 



The canonical value of Vsys for M51 is 472 ± 3 km s~^ ( iTullyl Il974l . and references 
therein). Table [2] lists TuUy's value for Vsys as well as the results from applying the three 
techniques described above. Our different methods give a mean systematic velocity of 471.7 
± 0.3 km s^^. Henceforth, we will fix Vsys to be 472 km s^^ (LSR, corresponding to a 
heliocentric velocity of 464 km s~^) in fitting the velocity field to estimate vr{iP) and V0{ip). 



4.3. Position Angle 

To investigate the effect of errors in the assumed galaxy position angle on the fitted 
values of vr and fe, we first use the aforementioned constant velocity model. Note that the 
position angle is required to deproject the galaxy image, as well as in equation ([2]). Using 
incorrect position angles, but correct model values of systematic velocity and inclination, 
yields a greater effect on vr than on ve. This is because an error ^6 ma <^ 1 in Oma 
results in an error —vg sin AO ma in the fitted vr, whereas the corresponding error in vg is 
~ +VRsm AOma- Since vg is large compared to vr, the shift in vr is larger than the shift in 
Vg. Thus , a ±10° error in position angle in the "constant velocity" model which has vg = 
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240 km produces approximately a =f40 km shift in the fitted radial velocity. Position 
angle errors also produce small perturbations in both velocity components. We conclude 
that using an accurate value of Oma is very important to obtain accurate vq and especially 
Vr fits. 

Unfortunately, the position angle of the major axis of M51 is particularly difficult to 
determine. The strong spiral arms and the tidal interaction with NGC 5195 distort the stellar 
disk, making it effectively impossible to determine a position angle from the orientation of 
the isophotes. Thus, it is necessary to go beyond morphology in determining 6 ma- We 
therefore re visit the determination of the galaxy's position angle. We apply the method of 



TuUyi (Il974l ) to our velocity data from Ha and CO observations. We also study the effect of 
streaming motions on this method by using model galaxies with known position angles and 
streaming velocities. In addition, we also apply two alternate methods to derive the position 
angle. 



The widely used value for the position angle of M51, 170°, was determined by iTuUy 



(jl974j ) using kinematic information. Tully assumed that the observed velocity should reach 
its extreme value at the position angle of the galaxy major axis, Oma- To determine OmAi 
he averaged the observed velocities in wedges extending over 5° in azimuth, and then for 
each radius took the position angle of the wedge with the extreme velocity as the estimated 
major axis at that radius. Tully excluded radii at which he was not confident that the true 
major axis had measured velocities (e.g. the faint interarm regions near the major axis). 

Figure [10] shows the results of applying TuUy's position angle determination method to 
Ha and CO observations. For each annulus of radial extent 5", the position of the wedge 
with the extreme velocity is marked. Due to the lack of data in the outer regions of the CO 
observations, only the wedges in the inner 70" provide reliable measures of the position angle 
of the extreme velocity. Similar to Tully, we did not attempt to estimate the position angle 
of the extreme velocity at radii for which data are sparse in the range of plausible position 
angles of the major axis. From the location of these extreme velocity wedges, we found Oma 
to be 172° (from an error weighted average) from both CO and Ha observations. 

However, streaming motions can shift the velocities, resulting in the extreme velocity 
occurring at position angles not corresponding to the true major axis. Indeed, inspection 
of positions of the velocity extremum wedges overlaid on the intensity maps shows that the 
position angles of the wedges in the interarms are clearly shifted counter-clockwise from 
those in the arms. This is most evident in the Ha maps, for which emission is detected from 
almost everywhere in the disk. We further explore streaming effects on the Tully method 
using a model with known streaming motions, generated using one of our vr and ve fits to 
the M51 CO data, with Oma = 170°. (Since this test is designed simply to reveal the twists 
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in the apparent position angle due to streaming, the particular value assumed for the true 
position angle and the particular vji and vg fits used is not significant.) Figure [TT] shows the 
results of applying Tully's method to this streaming model galaxy. In any given annulus, 
the extreme velocity averaged in the 5° wedges occurs in the interarm regions. In Tully's 
analysis, however, only spiral arm regions (near the apparent major axis) were considered, 
due to observational limitations. Therefore, the major axis position angle he found is likely 
biased clockwise from the true major axis. As shown in Figure [TTl even if the interarm 
regions are considered, the position angles of the locations of the extreme velocities do not 
necessarily correspond to the major axis. Thus, such an extrema method can be biased due 
to the inherent streaming in M51, regardless of whether the arms or interarms are considered. 

We employ two alternate position angle determination methods that make use of the 
full observed velocity field. In the first method, we average the observed velocity at each 
position angle in a wedge for both the northern and southern sides of the galaxy; then we 
fit a cosine curve to these averaged velocities as a function of azimuthal angle. We will refer 
to this method as the "radial-averaged" method. This is most easily accomplished in the 
polar projection, where we can average along a column to perform the radial average. The 
radial-averaged velocity as a function of position angle is shown in Figure [T21 along with the 
corresponding cosine fitted curves. We assume that the galaxy major axis should be at the 
position angle of the extrema of such curves. The mean position angle determined from the 
Ha and CO fits is ~177°, larger than the position angle determined by the Tully method. 
Again, the position angle determined in such a way is sensitive to streaming. Earlier we 
showed that streaming tends to cause the position angle of extreme interarm velocities to be 
biased counter-clockwise from the true value. Since the interarms occupy a greater fraction 
of the galaxy compared to the arms, streaming will introduce a counter-clockwise bias to the 
apparent position angle of the major axis. The effect of streaming on this method is further 
discussed below, following a discussion of our second position angle determination method. 

Our second method to determine the position angle, the "azimuthal fit" method, is 
similar to the one described in the previous paragraph, but instead the cosine curve is fit to 
the observed velocities along a projected circle with constant (projected) radius. As in the 
previous method, the polar projection of the velocity field is useful; in this case we simply fit 
a cosine curve to the velocities along a row of constant projected radius. The results from 
applying this method are shown in Figure [131 Note that the position angle of the velocity 
extrema varies systematically as a function of radius; it is approximately 180° in the inner 
region 30" from the center, declining to 165° 120" from the center; this trend, including the 
rise near Log(i?) = 1.7, is also evident from simple inspection of Figures S] and O Averaging 
over the radius range displayed in Figure [131 we obtain the same position angle of ~177° 
as in the previous method. Again, the velocities in the interarms bias this determination of 
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position angle, for the same reasons stated in the previous paragraph. 

In order to understand the effect of streaming motions on the position angle of the 
major axis derived using the radial averaged and azimuthally fit methods, we apply these 
methods to streaming models with known position angles and streaming velocities. The 
model velocity fields are produced from our vr and vq profiles obtained by assuming fixed 
values of 6 ma] "we then apply the radial averaged and azimuthally fit methods to these model 
velocity fields. Both methods recover a 9ma of ~176°, for all models, similar to the actual 
M51 velocity field, even though the position angles assumed in generating the streaming 
models can be very different. This is because for different position angles, the streaming 
velocities were derived from fits designed to best match the observed velocities. So in fact all 
the streaming models give virtually identical observed velocity fields regardless of assumed 
position angle. 

However, if we recreate the models setting the radial component to be zero everywhere, 
then we correctly recover the assumed position angles. This is clear evidence that radial 
streaming affects methods to determine 9m A) not only near the minor axis, as recognized by 
Tully, but also elsewhere including even the major axis. 

In order to quantify the effect of the non-zero radial velocities on the apparent position 
angle, we apply the position angle determination methods to models with known constant 
radial velocities. In these artificial models tangential streaming velocities are assumed to 
vary with arm phase, but radial streaming is assumed constant. We found that for every 
±10 km s~^ in radial velocity, the derived position angle differs from the actual position 
angle by ± 3°. This degeneracy between the position angle and radial velocity renders it 
difficult to accurately identify the true position angle, or to map the radial velocity. In order 
to accurately determine the position angle, we need to know the radial streaming. But in 
our effort to map the radial and tangential velocities of M51, we need to know the position 
angle. Thus, as we carry out our investigation, we shall use a range of position angles in 
deriving the two dimensional velocity components of M51. 

4.4. Inclination 

Estimating the inclination based on the orientation of the isophotes is unreliable, as 
discussed in §4.3[ due to the strong perturbations from the spiral arms and the tidal interac- 
tion. In principle, the inclination can be determined from a fit to the velocity field, as we did 
to obtain the systematic velocity in §4.21 However, the fit inclination is not well determined 
by the available data, presumably due to the streaming. To test the sensitivity of the vr 
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and ve fits to the inclination, i, in equation ([2]), we assume incorrect values of the inclination 
which differ from the true value by Ai for the constant velocity model (where vr and 
are constant). We then fit for vr and vq. As expected, an error Az results in errors in the 
fit value of vr and vq with magnitudes oc sinAz, along with small perturbations about this 
offset. 

Since we find that error introduced in the velocity components due to an incorrect 
inclination can be large (although just a sim ple scaling), we soug ht other constraints on 
the inclination. In particular, the TuUy-Fisher (ITully fc Fishen 119771 ) relation can be used to 
estimate the inclination. The well known Tully-Fisher relation is a correlation between galaxy 
luminosity and maximum rotation speed. The inclination can be estimated by comparing 
the rotational velocity predicted by the Tully-Fisher relation with the observed velocity of 
the fiat part of the rotation curve. W e use the baryonic for m of the Tully-Fisher relation 
discussed bv iMcGaughl fcoosh (see also lMcGaugh et al.l[200o[ ): 



M. = 50 



(3) 



where is the baryonic mass (in Mq), and Vc is the circular rotational velocity (in km 
s~^). Since the dispersion in the Tully-Fisher relation, L oc V^, is relatively small, given the 
luminosity the uncertainty in Vc is small. 

In order to determine the baryonic mass M;,, we require the stellar mass M^, which is 
related to the 5-band luminosity Lb and the B-band mass to light ratio (M/Lb), 



M, = Lb- {M/Lb). 



(4) 



We use the c orrelation of the galaxy co lor with M/ L discussed by iBell fc de Jond (120011 ) 
applicable to the lCharlot &: Bruzuall (Il99ll ) population synthesis models, to obtain {M/Lb): 



{M/Lb) = 10t-"-^^+i-5^^^-^^J. (5) 
The RC3 catalog Jde Vaucouleurs et al.lEggih gives {B-V) = 0.53 for M51, so {M/Lb)=IM. 



The last quantity required to determine M^ is the luminosity Lb, which can be de- 
rived if we know the distance. Two independent studies have given similar M51 distance 
estim ates: observation of p lanetary nebulae gives a distance modulus of m — M=29.62 ± 
0.15 ( Feldmeier et al.lll997l ). and a study of surface brightness fiuctuations in the companion 



NGC 5195 gives m — M=29.59 ± 0.15 (I Jensen et al.lll996l ). We thus employ a distance 
modulus of m — M= 29.6, corresponding to a distance of 8.4 ± 0.6 Mpc. Using the RC3 
catalog value of B = 8.67, corrected for extinction. 



^Q-0.4(8.67-29.6-5.48)^^ = 3.66 X lO^^L^. 



^0. 



(6) 
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Using this value and the {M/Lb) value of 1.54 (equation [5]) in equation (jl]), we obtain 



M, = 5.64 X 10^°Mc 



0- 



(7) 



We can now apply the iMcGaughl (120051 ) relation in equation ([3]) to obtain the circular 
velocity: 

K = [(M, + M,„,)/50]i/^ = 188 kms-\ (8) 

where Mgas is the total gas mass; in the case of M51 the gas is predominantly molecular. From 
our CO observatior is, we compute M„ as = 5.4 x lO^M©, using an X-factor of of 2 x 10^° cm~^ 
[K km s~^]~^ (e.g. IStrong et al.lll988l ). Due to the small value of the exponent in equation 
[SI errors in the mass, due to variations in the X-factor, for example, will not significantly 
affect the resulting rotational velocity. 



The observed velocity is related to the circular velocity by 

VcMbs = Vc smi. 



(9) 



Adopting the center, 6 ma, and Vgys described in this section, we apply a tilted- ring analysis 
to determine the flat part of the rotation curve. We obtain an observed circular velocity 
between 70 and 80 km s~^, implying 



22° < I < 25° 

Therefore, for our subsequent fits, we adopt an inclination of 24°0 



(10) 



4.5. Summary: System Parameter Values 

In summary, we have shown that the fit values of vr and vq are sensitive to the assumed 
values for the fixed parameters in equation ([2]), V^^^, Oma, and i. Uncertainties in the assumed 
position of the dynamical center are too small to significantly affect the derived streaming 
velocities. We have used three different methods to determine Vgys, which resulted in a value 
similar to the Vsys found by Tully. We have found it to be extremely difficult to constrain the 
value of the position angle of the major axis using the velocity field, due to the significant 
streaming that shifts the position angle of the extreme velocities. As a result, in fitting for vr 
and fe, we allow for a range of plausible position angles. Lastly, we will adopt an inclination 



-'^ Employing the standard TuUy-Fisher relat ion instead, we obtain a mean inclination of ^^23°, using the 
slope and zero-point fits from lVerheiienI (|200ll ). 



- 16 - 



of 24 °, which is determined by using the baryonic Tully-Fisher relation (IMcGaugh et al. 



2OOOI ) between the baryonic mass and rotational velocity. To estimate vg and vr, we thus 



use the center position and systematic velocity listed in Table [H but use a range of position 
angles and an inclination of 24°. 



5. Results: Velocity Profile Fits 

With our improved estimates of the global parameters, we apply the fitting algorithm to 
the observed velocity field in different annuli to determine the radial and tangential velocities 
vr and vo as a function of arm phase ip. We initially adopt a position angle 6ma of 170°. We 
address the issue of a varying 6m a in §6.2[ Figures [H] and [15] show the CO and Ha vr and 
ve fits in 6 overlapping annuli between a galactocentric radii of 21" and 105", and Figure [16] 
shows the corresponding overlapping annular regions. 

An initial inspection of the streaming profiles indicates that the velocity structure is 
rather complex. Models of density wave streaming qualitatively predict that as gas encoun- 
ters the arm, vr, which was positive (i.e. outward) in the interarm, becomes negative, and 
as the gas exits the arm again becomes positive. The azimuthal velocity vg is predicted to in- 
creas e rapidly as gas flows thro ugh the arm, and then decline more gradually in the interarm 
(e.g. [Roberts fc Stewartlll987l . and §31). First we concentrate on Arm 1 (the brighter arm, 
shown at '?/' = 0° and more fully at 360°). For vr there is a pronounced minimum close to the 
arm position, seen in both CO and Ha. There is a boost in vg through the arm, again seen in 
both CO and Ha. For this arm, the streaming is qualitatively as expected from steady-state 
spiral shock models. The velocities associated with Arm 2 (located aX ip ^ 200°) however 
do not agree with simple predictions. For vr, a clear minimum is only apparent in the outer 
annuli, and the boost in vg is weak or nonexistent. In the interarms, the structure appears 
somewhat more complex than the simple model expectation of a relatively constant or slowly 
rising vr and a slowly declining vg. We suggest that Arm 1 matches simple theory because 
its structure is simple, i.e. well described as a log spiral of constant phase. By contrast, for 
Arm 2 the CO distribution is not as well described by a single log spiral segment. Instead, 
it has several segments with different pitch angles and jumps in phase; thus the velocities 
associated with this arm is complex. 



One explanation for the differences in the two arms, as discussed by [Rix fc Riekd (Il993[ ). 



is that the spiral pattern in M51 is actually a superposition of a stro ng spiral mo d e wit h 



a m=2 Fourier component with weaker m=l and m=3 components. [Henry et al.l (|2003[ ). 
using the spatial distribution of CO emission obtained from the BIMA CO map, found such a 
scenario to be feasible by explaining the bright arm as the result of constructive interference 
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between the m=2 and m=3 components, and the weak arm the result of a destructive 
interference between the two components. There is evidence for interarm structure possibly 
supporting such a multiple density wave component description of the spiral arms, which 
would be expected to manifest itself in the kinematics. Recent Spitzer observations of M51 
clearly show spiral structure between the main CO arms; the Spitzer image and interarm 
features are discussed in the next section. 

There are also clear differences between different annuli. For example, in the 36" - 61" 
annuli, the Vr increase downstream from the weaker arm is much more pronounced than in 
the 21" - 36" annuli. In addition, there is a vr decrease to as low as < —50 km s~^ in the 
arms of the outer regions, which perhaps can be attributed to an incorrect choice of a fixed 
position angle for the disk (see §4.31) . There are also clear differences in the ve gradients 
between different annuli. 

By and large, similar velocity structure is apparent in both CO and Ha. For example, 
in the 47" to 80" annulus, the gradual rise in vr from -50 km s~^ at ip=l'^Q° to 70 km s~^ at 
■?/'=300° is shown in both tracers. Further, there is a strong vr peak at il)=12Q° in the 27" to 
47" annulus in both CO and Ha; however, such pronounced local extrema in the interarms 
are not expected in the theory for a single spiral mode. In general, the overall amplitude of 
the streaming and the location of most features coincide, and regions in which the velocity 
structure is somewhat different tend to be interarm regions where little CO is detected. 

Such similarities are not unexpected due to the dynamical coupling between the different 
components. CO, which traces the molecular component, is dynamically cold, with a velocity 
dispersion of only 4-8 km s~^. Thus, the molecular component of the disk reacts strongly to 
any perturbation, as evident in the strong vr and ve gradients associated with the spiral arms. 
The spiral arms compress the gas, triggering star formation. The newly formed hot O and B 
stars subsequently ionize the surrounding gas, resulting in Ha emission. Due to the fact that 
much of the Ha emission comes from gas near the region of birth, observed Ha velocities 
will be similar to observed CO velocities. There may also be a diffuse ionized medium not 
closely associated with the O and B stars, and this medium is likely not dynamically coupled 
with the molecular gas. However, as can be seen in Figures H] and 0, the brightest regions of 
Ha emission occur just downstream from the molecular spiral arms. Further, the generally 
good agreement between the velocity measurements from CO and Ha observations suggest 
that most of the ionized emission originates in gas associated with star forming regions. This 
similarity in velocity structure derived from independent observations also gives confidence 
that the fitted velocities are reliable and that the deviations from simple theory, including 
interarm features, are real. 

The profiles in Figures [H] and [15] qualitatively agree with previous studies of streaming 
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in M 51 involving ID cuts along the major and minor axes (e.g. iRandl Il993l : lAalto et al. 



19991 ). The gradients of the velo city profiles throu gh the arms in different annuli is in 



accordance with the conclusion of Aalto et al. Jl999h supporting the presence of shocks in 
the arms from a qualitati ve comparison of velocities along ID cuts to streaming models of 
Roberts fc StewartI (119871 ). In §6] we analyze the feasibility of a steady or quasi-steady spiral 



pattern in M51, which has been a working hypothesis for many analyses of the spiral arms 
of this galaxy. 



5.1. Interarm Structure 



In estimating the radial and tangential velocity components, we fit observed velocities 
along log-spiral segments. The slope of the slit is determined by the slope of the main CO 
arms on the logarithmic polar projection, i.e. the pitch angle of the arms. Though the slope 
of the CO arms, or at least the bright arm, is well defined, that slope may not be appropriate 
for the interarms. In other words, velocity may not be constant along the interarm log-spiral 
segments congruent to the main CO arms. 



The recent Spitzer 8 /im image of M51 (jCalzetti et al.l l2005l : iKennicutt et al.ll2003l ). 
shows clear interarm features not seen in the CO map due to the lower resolution of the CO 
observations. Many of these featu res are spurs (or feat hers) which have been found to be 
ubiquitous in grand design spirals (ILa Vigne et al.ll2006l ). These interarm features will also 
cause kinematic perturbations. In fact, close inspection reveals that interarm perturbations 
in the velocity field of M51 coincide with strong interarm features apparent in the 8 fim 
image. Since the features have different pitch angles from the main CO arms, we are likely 
smearing out these finer interarm velocity perturbations. As a result, the interarm velocity 
profiles we have derived do not reveal the details of the velocity perturbations associated 
with interarm substructure; in a detailed study of the 2D velocity field between the main 
arms, the interarm structure would need to be considered. 

In the next section, we use the vr and vg fits to assess the feasibility of the hypothesis 
of a quasi-steady pattern. Again, our fitting method is designed to reveal streaming solely 
associated with the spiral arms, and does not capture sm aller scale perturbations, such as 



those associated with interarm fe atures. Both observations dElmeg: 



reen 



1980 



La Vigne et al. 



20061 ) and numerical simulations (IKim fc Ostrikeiil2002l : IShetty fc Ostrikerll2006l ) have shown 
that spurs and feathers are associated with star formation, indicating that these features 
are not long lasting. The modal theory, hypothesizing quasi-stationary grand design spiral 



struc ture, acknowledges that such smaller scale features can be transient (e.g. iBertin &: Lin 
19961 ). In this study, smearing out the interarm perturbations likely does not affect the 
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overall conclusions we draw from the fitted velocity profiles. 



6. Tests of Conservation Laws 

6.1. Conservation of Vortensity 

For a flattened system, the conservation of mass and angular momentum can be com- 
bined to yield 

m [ E ) + ■ ^ [ E ) = °' ^'^^ 

where S is the surface density, and VmerUai is the velocity in the inertial frame. Equation 
(fTTj) states that the vorticity per unit surface density, known as vortensity, is conserved along 
streamlines. For steady systems, the conservation of vortensity can be simplified: 

^ ■ ( 1? + - 4^ ) = constant, (12) 



J: \R dR R 89 ^ 

because the temporal term in equation fllip vanishes. Even if the flow is not steady, portions 
of the galaxy that originated in a region of constant vortensity will still satisfy equation [T2j 

In order to test whether equation f|T2l) is satisfied for the gas in M51, we need the surface 
density E, which we can estimate using the observed CO brightness to derive the correspond- 
ing H2 column densityEI Most studies suggest that the relationship between CO and H2 is 
reasonably linear, though the conversion factor, known as the X-factor, is controversial. In 
our analysis, we will assume that CO is in deed a linear molec ular tracer, and employ an 



X-factor of 2 X 10^° cm'^ [K km s'^]"^ (e.g. IStrong et al.lll988r ). 



We first test the vortensity condition from the velocity profiles derived in the 47"-80" 
annulus, (see Figure [H] - [16]) . We choose this particular annulus because ve variations are 
relatively smooth in both arm and interarm regions, and are likely due primarily to spiral 
streaming. This annulus clearly shows the characteristic vg boost in the arm, and the more 
gradual interarm decrease in vg. We consider the fits derived from this annulus assuming a 
9m A = 170° (see Fig. [12]). As shown in §4.31 changes in 9m a affect vg only modestly. 

For the first term in equation (fT2|) . vg/R, we use the mean value of the tangential 
velocities fit in the given region. To measure the radial gradient of the tangential velocity. 



^As we will describe in ij6.2.21 HI can be neglected since the gas in M51 is mostly molecular in the region 
studied. 
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which appears in the second term, we use 

_d_ 
dR 



Rtanip dip 



" (13) 



where ip is the pitch angle of the spiral arms, and we assume negligible variation parallel 
to the arm. We adopt a pitch angle of 21.1° (which is also the slope of the "slit"). We fit 
straight lines to the velocity profiles in order to approximate the last two terms in equation 
f|T2|) . For the surface densities, we use the peak value for the arm, and for the interarm we 
use the value of S at phase separated by 90° from the arm. Again, we are assuming that CO 
directly traces the molecular abundance. Table [3] shows the vortensity values for the 47"-80" 
annulus, including the values of each of the terms in equation (fT2l) . 

Table[3]shows that, within the errors, the arms and Interarm 2 have consistent vortensity 
values. The value for Interarm 1, however, is lower than in the other regions. The lower 
value for Interarm 1 can be inferred directly from the profile itself (Fig. [T^ . The tangential 
velocities in both arms are clearly rising, and Vg in Interarm 2 (downstream from Arm 2 
at '?/' = 200°) is predominantly decreasing, suggesting that the spiral arms have the most 
significant influence on the velocities in these regions. In Interarm 1 (downstream from Arm 
1 at = 0°), however, there is more structure to the velocities, suggesting that there are 
other sources of perturbations in addition to the spiral arms. We find that for most of the 
velocity profiles in Figure [H] and [TS] the vortensity values are consistent between the two 
arms, within the errors. However, in the interarms, the vortensity values differ. We find 
varying vortensity values in all interarm regions except for Interarm 2 indicated in Table [31 
Overall, the agreement between vortensity in arm regions in each annulus indicates either 
that a steady state depiction of the vortensity is valid and there is very little radial migration 
of gas, or else that in a given annulus much of the gas originated in a region of constant 
vortensity and has been conserved along streamlines as gas in a whole annulus flows inward 
or outward. 



6.2. Conservation of Mass 



6.2.1. Flux Weighted Average vr 



In the QSSS scenario, the spiral pattern — as defined by its amplitude, phase, and 
rotation rate — would not change significantl y over the cour s e of a few revolutio ns in a 
frame rotating along with the spiral pattern (ILindbladl Il963l : iBertin &: LinI Il996l ). Such 
a framework suggests that on average any accretion of material into the arms should be 
balanced by the same amount of material exiting downstream. This condition corresponds 
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to conservation of mass for a steady state system; if this condition holds it should be apparent 
in the variation of observed velocities with spiral arm phase. 

As can be seen in Figures [T3] and [TSl the spiral arms clearly perturb vq with deviations 
of > 100 km s^^ but vg always remains positive, indicating that the orbital flow is in one 
direction only. However, the radial velocities do change sign, indicating both inflow and 
outflow. If the spiral arms are indeed a quasi-stationary pattern, then large amounts of 
matter should not be undergoing net inflow or outflow, i.e. the mass-weighted average radial 
velocity cannot be too large. A large or spatially strongly variable mass-weighted average 
vr would imply a very dynamic system. In particular, if the sign of this quantity changes, 
then there would be a buildup or depletion of mass in one or more radial locations. 

As discussed in §4.31 the fit for vr is very sensitive to the assumed value of the position 
angle. Thus, in investigating the mass-weighted average radial velocity, we consider a range 
of position angles. Figure [17] shows the vr fits for Ha in annuli with radii of 47" - 80" 
and 61" - 105" for three different position angles, 170°, 175°, and 180°. The CO fits are 
similar, but noisier and have larger error bars (see Figs. [T^ - ITH) . A striking aspect of the 
fits in Figure [T7] is the large magnitude of inflow in the arms for all three position angles; 
the radial velocity drops to as low as —75 km s~^, suggesting significant inflow for gas in 
the spiral arms. In the upstream regions vr is positive, approaching 70 km s~^ for some 
parameter choices. For a region farther in, in the 27" - 47" annulus shown in Figures [H] and 
[T^ the fitted radial velocity (assuming 6ma = 170°) reaches values greater than 100 km s~^, 
indicating tremendous outflow in the inner regions; assuming a position angle of 180° for 
this annulus only reduces the peak velocity from RilOO km s~^ to ~75 km s~^. 

Figure [18] shows the flux- weighted average radial velocity, (f/j), for different position 
angles in the different annuli used in the fitting process. With the canonical position angle 
of 170°, there is significant outflow in the inner regions of M51. On the other hand, a 
position angle of 180° seems appropriate for the innermost regions of M51, since this yields 
a lower value of {vr). However, with such a position angle we find significant inflow in the 
outer region. If we adopt an intermediate position angle of 175°, there is outflow in the inner 
regions and inflow in the outer regions. For a 9ma of 175°, (vr) = for the 36" - 61" annulus, 
with mean radius (R) = 1.98 kpc, while adjacent annuli have (vr) = 10 and —20 km s~^, 
for (R) = 1.5 kpc and 2.6 kpc respectively (1" = 40.7 pc at a distance of 8.4 Mpc). If this 
were true, then the gas would all collect near R ^ 2 kpc in less than one orbital time-scale 
(~200 Myr), which is not consistent with a steady state. 

This analysis leads us to conclude that if the spiral pattern is long-lived, the large 
variations in the radial velocity shown in Figure [18] suggests that the position angle must 
vary with radius, indicating a disk which is not coplanar. This trend suggesting a larger 
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position angle in the inner regions and a smaller position angle in the outer regions is also 
in accordance with the position angle tests described in §4.31 (see Figure [T3|) . 

We schematically show a disk with a varying position in Figure [T^ the position angles of 
the ellipses are arranged as indicated by Figure [131 As discussed, one effect of a variation of 
position angle is a disk that is not coplanar. The inclination in this schematic is exaggerated; 
the observed morphology, including the apparent spiral structure, depends on the viewing 
angle, among other factors. 



6.2.2. Continuity and Spiral Pattern Speed 

In this section we explore the plausibility of QSSS using the gas continuity equation. 
The continuity equation for gas flow in a two-dimensional system is 

— + V ■ (Ev) = 0; (14) 

this holds in any frame, e.g. whether the velocity is measured in an inertial frame or one 
rotating at a constant pattern speed. The first term, dH/dt, represents the temporal growth 
or decay of the surface density E at any given radius R and azimuthal angle 6 in the plane 
of the galaxy, where those coordinates are with respect to the frame in which the velocity is 
being measured. 

If the flow is in a steady state, then the temporal term vanishes, leaving only the mass 
flux term Sv. If the gas is responding primarily to a single dominant spiral perturbation, as 
would be required for a fixed spiral pattern, and when v is measured in the frame rotating 
at the pattern angular velocity Qp, 

V ■ [n^^nerUal ' ^pRO)] = 0. (15) 

Thus, for an exact steady state the mass flux must be constant (in the frame rotating with the 
same angular velocity as the spiral mode). For a quasi-steady state, the temporal variations 
in E will only be small, and thus variations in mass flux would also be small. This condition 
can be further simplified using a reference frame aligned locally with the spiral arms. Figure 
[20] shows this reference frame; the x and y coordinates are the directions perpendicular 
and parallel to the local spiral arm, respectively. The transformation between cylindrical 
coordinates and this arm frame is achieved using 

X = cos ipR + sin ip9, (16) 



— sin ipR + cos ip9, 



(17) 
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where ip is the pitch angle of the arms. 

The velocities in the arm frame are given by 

Vx = VRCOsip + {vg — QpR) sinip, (18) 

and 

Vy = —vr smip + {vq — VLpR) cosip. (19) 

From the maps shown in Figures [T] and HJ it is apparent that the intensity and ve- 
locity vary significantly more across the arms, in the x-direction, than along them, in the 
y-direction. Thus, the variation in the product of Hvy along y is much smaller than the 
variation in the product of Ev^. along x, reducing equation ([T^ to 

I]Vx ~ constant. (20) 

Namely, for a steady pattern, as the gas decelerates (in the x direction, perpendicular to the 
arm), mass accumulates and the surface density increases; as the gas velocity increases, the 
surface density decreases. 

One difficulty in testing whether equation fl20l) holds is that neither 6ma nor Qp is well 
constrained. Errors in 6m a yield errors in the fitted value of vg of ~ vr sin A6ma and in the 
fitted value of vr of ~ —vg sin A9ma- If Afip is the error in the pattern speed, then the fitted 
value of Vx will be approximately given by Vx + {— vg sin A9 ma cosip + VRsin A9 ma sin ip — 
AQpR sin ip) . Since the vr term has two factors of the sin of small angles, that term will be 
much smaller compared with the other two terms. The true value of Vx will therefore differ 
from the fitted value by Cx ~ vg sin AO ma cos ip + AVLpR sin ip. 

In order to assess whether steady state continuity as expressed by equation fl20l) holds 
in the case of M51, we therefore consider the quantity 

YiVx = T,{vRCOsip + {vg — QpR) sinip + C), (21) 

where vr and vg are fitted values and C = {Cx), i.e. the (unknown) azimuthally-averaged 
correction due to the errors in 9m A and flp. We apply equation (pT]) by solving for the value 
of C using the values of Vr, Vg, and E in the two arm segments of an annulus: 

^ _ [vR,arml COS ip + {vg^arml " ^pR) sin ip]Earml — [vR,arm2 COS ip + {vg 

,arm2 ^p-R) siu 2p] Efjrm2 

" y Tv 

(22) 

We then test whether the value of C obtained using equation fl22|) also satisfies equation fl2Tl) 
in the interarm regions. If equation fl2Tl) is satisfied for both interarm and arm regions, it 
would suggest an approximate steady state. 
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We again focus on the 47"-80" annulus, where the vr (and vg ) are relatively "smooth." 
The vr profiles for this annulus is shown in Figure [T71 assuming three different values of the 
position angle. Table H] shows the relevant values associated with equation fl?Il) for the arms 
and interarm r egion s. We employ the pattern speed of 38 ± 7 km s~^ kpc"'^, calculated by 



Zimmer et al.l (120041 ) by applying the Tremaine-Weinburg method CO observations of M51. 
Corotation corresponding to this pattern speed is marked on Figure [161 After solving for C 
using quantities from the arms, it is clear that the flow in the interarm region is not consistent 
with a steady state description. A position angle of 170° produces large negative mass flux 
in the arms, and positive flux in one interarm region. Even variations in the X-factor cannot 
resolve this discrepancy. Assuming larger values of the position angle still produces mass 
fluxes with vastly different magnitudes, and even different signs. Increasing the error in C 
up to an order of magnitude still cannot result in consistent mass fluxes between the arm 
and interarm. This suggests that any reasonable changes to the values of 9m a or Qp will still 
result in varying mass fluxes. We have checked the mass flux in other annuli using the same 
method as for the 47"-80" annulus, as well as in other localized regions not presented here, 
and found similar discrepancies in the mass flux. 

In our analysis of continuity so far, we have not taken into account the contribution 
from the atomic component of t he disk. In fact, in mo st galaxies the majority of the gas 



exists in the form of HI. In M51. iTilanus fc Allen! (119891 ) showed that the downstream offset 



of HI relative to the dust lanes is likely due to dissociation of molecular gas by recently 
formed massive stars. However, the inner disk of M51 has an unusually large fraction of 
molecular gas, so even at peaks of the HI photodissociation arms, the contribution of at omic 



gas to the total gas surface density is negligible. Using the HI maps of iRots et al.l (jl990l ). we 
find that the atomic column density N(HI) is significantly less than the molecular column 
density N(H2) in the vast majority of locations in the inner disk (21" < R < 105"); N(HI) 
exceeds N(H2) in only ~ 7% of the inner disk. The mean value of N(H2)/N(HI) throughout 
the inner region is ~10. Though we used a constant X-fact or to obtain the molec ular surface 



density, moderate variations in the X-factor (for M51, see iNakai &: Kundll995l ) will not be 
sufficient to account for the discrepancy. Nevertheless, no change in the X-factor, or in the 
contributions of the molecular or atomic matter to the total mass, can account for the change 
in sign of the mass flux; the varying sign of Ef^, can only be due to a sign change in v^, not 
E. 

Our conclusion, after analyzing the mass flux, is that the kinematics are not consistent 
with a quasi-steady spiral pattern in a flat disk. We find that no single pattern speed can 
satisfy quasi-steady state continuity, suggesting that the QSSS hypothesis is not applicable 
to M51. It is essentially the tremendous variations of the radial velocity within a given 
annulus — amounting to ~ 100 km s~^ — that lead the QSSS hypothesis into difficulty. 
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One explanation for the transient nature of the spiral arms in M51, perhaps due to the 
interaction with its companion, is a spiral perturbation with a constant pattern speed, but 
with time-varying amplitude. Or, there may be multiple modes at work in the disk of M51, 



which may be construed as a mode with a radially varying pattern speed (e . g. iMerrifield et al. 



20061). Multiple patte r ns spe eds in M51 have been previously sugg ested bv lVogel et aP Jl993h 



and lElmegreen et al.l (119891 ). However, the extreme variations in the radial velocity cannot 
be explained by multiple patterns alone. Possible causes for the large observed vr gradients 
are large out-of-plane motions or a variation in inclination; since the inclination of M51 is 
small, a variation in i due to a warped or twisted disk will produce large variations in the 
observed velocity due to projection effects. 



6.3. Discussion 

We have shown that the density and velocity structure in M51 does not support a quasi- 
steady state depiction for the spiral pattern, using measurements of the mass flux. Further 
evidence that the observed structure is inconsistent with steady state can be obtained by 
adopting the fitted 2D velocity field, and demonstrating that the density structure is then 
non-steady . We have c arried out this exercise using a modified version of the NEMO task 



FLOWCODE ( lTeubenlll995l ). In this exercise, a disk is populated with gas tracer particles using 
the intensity profiles averaged along spiral segments, reproducing the spiral density pattern 
of M51. Each location in the disk has an associated vr and vg, given by the fitted velocity 
profiles (e.g. Figures and [T5|) . and an assumed value of the pattern speed. The motion of 
the particles is then integrated using FLOWCODE: after a suitably small timestep, the particles 
take on new velocities depending on their location in the disk. In essence, this simulation is a 
purely kinematic test to determine whether the steady state continuity equation (eqn. [T^ ) 
is satisfied or not, using the density and fitted velocity profiles of M51 (Figs. [T^lfTSj) . We 
find that the input spiral pattern vanishes in less than one orbital time scale (~200 Myr), 
regardless of what values of the position angle and pattern speed we assume. 

The precise nature of the velocities is one of a number of issues that need to be considered 
in further studying the the global spiral pattern in M51. For example, our result suggests 
the role of a warp certainly needs to be taken into account. There are strong indications 
that the outer disk of M51 is warped; our finding suggests that the disk is not coplanar 
even further inward. The non-coplanar attribute may be the result of the tidal interaction 
between M51 and its companion. 

The possible warp and/or twist in the disk of M51 would of course affect the projected 
velocities, and would present itself as gradients in the velocity components, as discussed 
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in the previous section. If this were indeed the case, then the single or multiple in-plane 
modes would have to be in phase with the vertical mode in order to sustain a spiral pattern. 
The inherent uncertainty in deriving three velocity components from the single observed 
component leads to difficulty in estimating and analyzing both the vertical and in-plane 
modes. 



7. Summciry 

We have analyzed the velocity field of M51, using CO and Ha observations, to investigate 
the nature of the spiral structure. We summarize the main results here: 

1) The velocity field is quite complex. Observed velocities show significant azimuthal 
streaming associated with the spiral arms, as well as strong gradients in the radial velocities. 

2) The aberrations in the velocity field strongly suggest that the disk is not coplanar, 
perhaps as far in as 20" (~800 pc) from the center. 

3) We obtain fitted radial and tangential velocity profiles by assuming that velocities 
in any annulus vary only with arm phase. Strong gradients in the radial and tangential 
velocities are found in the profile fits. In general, the shape of both the vr and ve profiles are 
in qualitative agreement with theory of nonlinear density waves, and support the presence 
of shocks. 

4) In detail, the velocity profiles from different radial regions of M51 differ significantly. 
In addition, velocity profiles associated with the two arms also show differences in a given 
annulus. For the arm that is well described by a logarithmic spiral (bright arm), the as- 
sociated velocities are in good agreement with simple theoretical spiral shock profiles. For 
the other arm, which is not as well described by a logarithmic spiral, the velocities are more 
complex. 

5) The velocity profile fits from CO and Ha emission are rather similar, suggesting that 
most of the Ha emission originates from gas associated with star forming regions. 

6) When we assume a single value for the position angle of the major axis of M51 and 
inclination, we find that large amounts of material flows toward an annulus of intermediate 
radius, due to the large gradients and change of sign in the flux weighted average radial 
velocity. As a result, cither the position angle of the major axis or the inclination must vary 
with radius, suggesting that the disk of M51 is warped and twisted. 

7) We analyze conservation of vortensity, using the radial and tangential velocity proflle 
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fits. We find that vortensity is fairly consistent within a given annulus, indicating that the 
gas there all originated in a region of uniform vortensity. 

8) Using the equation of continuity, we find that the density and fitted velocity profiles 
are inconsistent with quasi-steady state mass conservation in any frame rotating at a constant 
angular speed, at least for a planar system. Variations in the pattern speed, position angle, 
and X-factor alone cannot account for the differences in the mass flux, suggesting that spiral 
arms are quite dynamic, and possibly that out-of- plane motions are significant. 

We are grateful to S. McGaugh for his help in using the baryonic TuUy-Fisher relation to 
estimate the inclination of M51. We thank W.T. Kim for useful discussions, and G. Gomez 
for providing velocity profiles from his hydrodynamic models. We also thank the anonymous 
referee for raising many interesting questions and useful suggestions. This research was 
supported in part by grants AST 02-28974 and AST 05-07315 from the National Science 
Foundation. 
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Table 1. Initially Adopted Parameters for M51 



Parameter Value Reference 



Center RA (a) (J2000) 13^29^52^71 

Center DEC (S) (J2000) 47°11'42'.'80 

Systematic Velocity (Vsys) 472 (LSR) 

Position Angle of Major Axis {9ma) 170° 

Inclination (i) 20° 



Table 2. Estimation of the Systematic Velocity of M51 



Method 


Vsys (km s-i)^ 


Error (km s ^) 


TuUy (1974) 


472 


3 


Freeing Vgys in fitting 


470.6 


0.4 


for Vr{iP) and Vg{ip) 






Tilted Rings Analysis 


471.4 


0.5 


Rotation Curve Fitting 


473.3 


0.5 


Weighted Mean 


471.7 


0.3 



Velocity in LSR frame 



Hagiwara et al. (2001) 

TuUv (1974) 
)) 

5! 
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Table 3. Vortensity in the 47" - 80" Annulus 



Region 






(dVg\h 

\ dRj 


( 1 dvii\h 

\R de J 


Vortensity Value'^ 


Arm 1 {ip^ 360°) 


244 


70 ± 1 


103 ± 2 


-30 ± 3 


0.8 ± 0.2 


Interarm 1 [ip ^ 90°) 


16 


67 ± 1 


-52 ± 3 


6 ± 1 


0.5 ± 0.2 


Arm 2 {tjj^ 190°) 


128 


63 ± 0.5 


59 ± 3 


-18 ± 3 


1 ± 0.2 


Interarm 2 {ip ^ 275°) 


19 


59 ± 1 


-31 ± 2 


11 ± 1 


0.9 ± 0.2 



^[M© pc-2]; error of ~20% 
^[km kpc^i (Mq pc"2)-i] 




Fig. 1. — CO (1-0) velocity-integrated intensity (bottom) and velocity (top) maps of M51. Velocity contours 
increment by 10 km s""'^, between 360 and 560 km s^^. Overlaid lines are logarithmic spirals with a pitch 
angle of 21.1°, separated by 180°. 
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Table 4. Mass Flux in the 47" - 80" Annulus 



Region 




vr cos ip^ 


(v0 — Qr,it!) sin ij}^ 






Oma = 170°: 












Arm 1 {tp ^ 360°) 


244 


-24 


41 


-58 ± 30 


-10113 


Arm 2 (ijj!^ 190°) 


128 


-53 


32 


-58 ± 30 


-10113 


Interarm 1 {ip 90°) 


16 


23 


18 




-259 


Interarm 2 (V' ~ 275°) 


19 


49 


21 




237 


Oma = 175°: 












Arm 1 ^ 360°) 


244 


-41 


41 


-32 ± 18 


-8103 


Arm 2 ^ 190°) 


128 


-63 


32 


-32 ± 18 


-8103 


Interarm 1 {ip ^ 90°) 


16 


12 


18 




-38 


Interarm 2 (i/j ^ 275°) 


19 


36 


21 




464 


Gma = 180°: 












Arm 1 360°) 


244 


-59 


41 


-8 ± 17 


-6345 


Arm 2 (ijj!^ 190°) 


128 


-71 


32 


-8 ± 17 


-6345 


Interarm 1 (V' 90°) 


16 


1 


18 




182 


Interarm 2 (?'• ?i 275°) 


19 


27 


21 




761 



^[Mq pc-2]; error of ~20% 
^[km s-^] 

■^[km s~^]; error largely due to errors in E and Vt. 
'^[Mq km s-^] 



-34- 




Fig. 2. — Hq velocity-integrated intensity (bottom) and velocity (top) map of M51. The overlaid spirals, 
as well as the velocity contours, are as described in Figure [TJ 
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< V > 

Azimuthal Angle 

Fig. 3. — Geometry depicting the spiral arm phase ip. The diagram on the left is the geometry in the 
plane of the galaxy. The phase, tp, represents the angular displacement between two locations with equal 
galactocentric radius Ro for two congruent spiral segments. The diagram on the right is the logarithmic 
polar projection of the geometry on the left, showing the corresponding spiral segments. 




yjH 7IK 'DO P 'iH 



Fig. 4. — Logarithmic polar projections of the CO intensity and velocity maps. Though the origin of the 
abscissa (azimuthal angle) is arbitrary, in this case it is aligned with due North. The direction of rotation is 
to the right (counter-clockwise as seen on the sky) . Also shown are the two logarithmic spiral lines positioned 
along the two spiral arms, which correspond to the lines overlaid on the maps of Figure [H 
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Fig. 5. — Logarithmic polar projections of the Ha intensity and velocity maps. Coordinate system and log 
spiral overlays are as in Figure [H 
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Arm Phase f (°) 

Fig. 6. — Gas profiles as a function arm phase ip for a model spiral galaxy (see text). Solid lines: density 
(top), Tjj^ (middle), and vg (bottom) profiles averaged at each arm phase, in an annulus extending from 8.38 
- 8.92 kpc. The vr and vg dashed lines in lower panels are obtained by fitting equation ([2]) to the "observed" 
velocities in the 8.38 - 8.92 kpc annulus. 
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Fig. 7. — CO (left) and Ha (right) vu and ve fits as a function of arm phase for an annulus with an inner 
radius of 21" and an outer radius of 36". The one sided 3a error-bars are also shown on the bottom of each 
panel. Dashed lines are the corresponding mean CO velocity integrated intensities, with the scale depicted 
on the right ordinate. Table [T] shows the fixed (canonical) parameters used in obtaining these fits. 
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480 




Arm Phase V (°) 



Fig. 8. — Result of fit to M51 CO velocity data in which Vgys, vji, and vg were allowed to vary. The radial 
range of the annulus is 14" - 136". The fact that the fit value of Vgys varies with phase tp shows that other 
parameters (e.g. i, 6ma) vary with radius within the annulus. The mean of the fits is 470.6 km s~^ (LSR), 
shown by the dashed line. 
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Arm Phase V (°) 

Fig. 9.— Fit Vs ys as a function of arm phase, similar to Fig. [3 but for a model galaxy with a warp. The 
systematic velocity adopted for the model is 400 km s^^ (dashed line), equal to the mean of the fits (solid 
line). 
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Fig. 10. — TuUy "wedge" method for estimating galaxy position angle. The extreme velocity for each 
5" annuhis, averaged in 5° wedges, is marked. The upper panels show the Ha velocity (left) and velocity- 
integrated intensity (right); lower panels show the same for CO. For CO, emission was too weak to apply 
the method at some radii, especially in the outer galaxy. 



Fig. 11. — Application of TuUy wedge method to a model with streaming. The velocity field is shown 
(left) along with the corresponding intensity map (right). The position angle assumed in the model is 170°, 
shown by the solid line. The 5° wedges with the extreme velocity for each annulus is marked. It can be 
seen that streaming shifts the estimated position angle from the true position angle. Note that the extreme 
velocities do not occur in the arm. 
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Fig. 12. — Mean observed velocity plotted vs azimuthal angle. All observed velocities were averaged over 
radius at each azimuth. Velocities are fit by cosine functions (solid line); extremum of the cosine curve 
indicates the best fit position angle of the major axis. 
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Fig. 13. — Fit position angle of the major axis 0ma as a function of galactocentric radius from Ha 
(upper) and CO (lower) velocity fields. The position angles were obtained by fitting a cosine function to the 
distribution of observed velocity vs azimuthal angle at each radius. 
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Fig. 14. — CO vr (left panels) and ve (right panels) fits as a function of arm phase tp 
different annuli (with radii labeled in the upper right of each panel). The thickness of the 
line shows a range of ±3a. Only vpt and ve fits with 3a < 20 km s~^ and < 60 km s~^, 
respectively, are shown. Dashed lines are the corresponding mean CO intensities, with the 
scale shown on the right ordinate. We assume a position angle of 170°, an inclination of 
24° and the center position and systematic velocity listed in Table [TJ Figure [16] shows the 
annular regions of M51 considered for these fits. 
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Fig. 15. — Ha VR and vg fits as a function of in different annuli, as in Figure [Ml 
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Fig. 16. — Deprojected CO map of M51 showing the overlapping annuh for which vr and vq are fitted 
as a function of arm phase (shown in Figures fT4l and [T5 | . The radii of (sohd) circles, from the inner to the 
outer, are: 21", 27", 36", 47", 61", 80", and 105". The annulus marked by dashed circles (4.2 kpc < i? < 
6.1 kpc) spans possible corotation radii corresponding to an adopted spiral pattern speed f2p = 38 ± 7 km 
s-i kpc-i (see ^6^21) . 



-48- 




Arm Phase {'4/) 

Fig. 17. — Ha VR fits as a function of arm phase for 3 different position angles Oma^ 170°, 175°, and 180°, 
for two annuli (47" - 80 " and 61" - 105"). We fix tlie inclination at 24°, and otlier parameters used in tlie 
fitting are shown in Table [T] The dashed line is the mean CO intensity along the arm for a Oma of 170°, 
which varies only slightly with Oma- The error bars are not shown because they are similar to those shown 
in Figure [T51 



- 49 - 



40 



20 



E 



A 

> 
V 



-20 - 



-40 




40 



20 



-20 



-40 



Fig. 18. — Mass -weighted average radial velocities (vr) in the different annuli, two of which are shown in 
Figure [T7l The abscissa indicates the mean radius (R), in arcsecs, of each annulus. The three panels show 
the mass- weighted average vn assuming three different values for the position angle 0m A- The error bars 
include both fitted errors in vn (see Fig. [14]) and an estimated error of 20% in E. 




Fig. 19. — Model disk showing the variation of the position angle with radius. The position angle profile 
is taken from Fig. [131 The inclination is exaggerated to show a more edge-on view. 
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Fig. 20. — Coordinate transformation geometry, from {R,6) galactocentric coordinates to the {x,y) spiral 
arm frame. 



